Non-perturbative Dynamical Decoupling Control and Threshold 
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Protection of quantum states in a noisy environment is essential for reliable quantum information 
processing. In this Letter, through a sequence of unoptimized periodical rectangular pulses we 
show that beyond perturbation, dynamical decoupling can be a powerful tool for the control of 
fidelity and entanglement of the open system in a non-Markovian bath. With this non-perturbative 
approach, we have found, for the first time, the threshold and pulse parameter regions in which the 
non-idealized pulses are effective in controlling quantum coherence. 
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Introduction. — The ultimate information processor in 
modern quantum technologies is a quantum computer 
While a full-fledged quantum computer is still im- 
practical at the present moment, research on quantum 
devices such as multiplc-qubit simulators 0, the quan- 
tum cryptography device 0), quantum imaging etc, 
has achieved momentous progress in recent decades. A 
simple, but vital quantum device called quantum mem- 
ory is a quantum system that can be used to store quan- 
tum states in a noisy environment for a certain period of 
time. Therefore, critical to quantum information process- 
ing is the ability to efficiently protect quantum memory 
from the corruptive influence of the environmental noises 
either classical [I[ or quantum 

A variety of theoretical strategies for combating the 
deleterious environmental noise have been proposed 0- 
15| . Among them, Bang-Bang control or dynamical de- 
coupling control of system-environment interactions by 
external fields is widely discussed due to its simplicity 
and wide applicability Although the theoreti- 

cal formalism for the dynamical decoupling is widely em- 
ployed in combating decoherence and protecting informa- 
tion leakage, most analyses and numerical simulations are 
based on idealized pulses and the Trotter product formula 
[ifj| . where it is assumed that these idealized pulses are 
so fast and strong that the system-bath Hamiltonian can 
effectively be turned off when the pulses are applied to 
the system. More specifically, if the control Hamiltonian 
is described by H c = Ja z , the total evolution operator 
in a short time 5 is U(S) = exp(— iH8 — iJSa z ) w —io z 
when J 5 = tt/2 and S goes to zero (The strength J goes 
to infinity). This amounts to taking the zero-order per- 
turbation for an effective Hamiltonian H c g — ir /2o~ z + 5H 
(setting h — 1) with the small perturbation parameter S. 
Without comparison with exact (numerical or analytical) 
solutions, it is not clear at all whether the zeroth-order 
perturbation or even perturbation theory is valid or not 
for any of models with realistic system-bath interactions. 



For instance, if a phase transition takes place when 5 is 
zero, then the omission of the system-bath interaction 
may become problematic because of the singularity at 
the critical point, no matter how small 6 is, exempli- 
fied by the well-known failure of perturbation theory in 
treating the BCS pairing model. Specifically, if we set 
one of modes in the BCS model as the system of interest 
and the others as bath modes [2lJ , the pairing correlation 
between this mode and others corresponds to the system- 
bath coupling and we discover a second-order quantum 
phase transition at 5 = 0. Therefore, it is highly de- 
sirable to have a non-perturbative theory to explore the 
non-perturbative regimes and provide the effective pa- 
rameter regions of the control pulse. 

In this Letter, we present a novel, non-perturbative 
dynamical control theory, based on the exact subsystem 
stochastic master equations obtained from combining the 
Feshbach projection technique [24] and the quantum 
state diffusion equation (QSD) [25-27]. Our approach 
employs the non-idealized pulses with finite widths A, 
periods t and finite pulse strengths $/A as characteris- 
tic parameters. Surprisingly, we have found the threshold 
and efficient parameter regions where the fidelity and en- 
tanglement dynamics can be controlled. We show that 
the characteristic quantity for quantum control is the ra- 
tio between the pulse period r and the width A. The 
other parameters, such as pulse strength, do not play an 
essential role in this scheme 23]. Our results show that 
the parameters for idealized fast-strong pulses control (A 
and r are effectively taken as zeros) occupy a very small 
portion of the parameter region in which efficient quan- 
tum control can be realized. In addition, we will show 
that the non-Markovian environment is crucial for the 
dynamical control of quantum coherence. 

Our non-perturbative approach can accommodate 
generic system-bath couplings, structured baths and var- 
ied numbers of environmental modes ranging from one 
mode, two modes to infinite number of modes. Our re- 
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suits provide higher flexibility in the experimental imple- 
mentation of the dynamical decoupling scheme since our 
scheme allows much larger parameter region. 

Stochastic master equation for subsystem — We start 
with the exact QSD equation 25|, 26, 28 -32| describing 
a quantum trajectory of an open quantum system ip t : 

d t ip t = [-iH sys + Lz* t - tfO(t, z*)]i/h = H eS (t)^ t , (1) 

where H sys and L are the system Hamiltonian and Lind- 
blad operator, respectively. z* is the Gaussian pro- 
cess @ describing the environmental influence, satisfying 
M[z*] = M[z*z*] = and M[z t z*] = a(t,s), where M 
denotes the statistical mean over the noise z t * and a(t, s) 
is an arbitrary environmental correlation function. 

Using PQ-partitioning technique, we can divide the 
n-dimensional wave-function ip t into two parts: an in- 
terested m-dimensional vector P(t) and the rest part of 
the vector, Q(t), which is an (n — m)-dimensional vec- 
tor. Accordingly, the stochastic state vector ipt and the 
stochastic effective Hamiltonian may be rewritten into a 
block form given by 
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Here the m x m matrix h and the (n — m) x [n — m) ma- 
trix D are the self-Hamiltonians living in the subspaces 
of P and Q, respectively. W and R are the coupling 
terms. For the selected m-dimensional subsystem-P, we 
can establish the following stochastic master equation: 



id t P(t) = h(t)P(t) + R(t)G(t, 0)Q(0) 



i I dsg(t,s)P(s) 
'o 



(3) 



where G(t, s) = %- |exp —i ds'D(s') j and g(t, s) — 

R(t)G(t,s)W(s). When m = 1, P is the amplitude of 
one basis vector, and Eq. ([3]) becomes a one-dimensional 
stochastic master equation As shown below, by 

choosing a proper basis and the order of the basis vectors, 
one can significantly simply the numerical simulation and 
analytical calculations. 

Models and Analysis — Here we consider two models: 
(i) a single-qubit system in dissipative bosonic environ- 
ment with H sys = -|cr z and L = cr_; (ii) a two-qubit sys- 
tem in a collective dissipative bosonic environment with 
H sys = §{o-^ + af) and L = er 4 +cr?, where E may be a 
time-dependent function. The former corresponds to the 
rather general individual error considered in quantum in- 
formation, and the latter to the collective error. 

For simplicity, the non-Markovian environmental noise 
throughout the paper is described by the Ornstein- 
Uhlcnbek type noise with the correlation function 
a(t,s) = Sre -7 '* -s ', where I/7 characterizes the memory 
time of the non-Markovian environment. The correlation 
function has a well-defined Markov limit when 7—5-00. 



For the single-qubit model, in the basis {|0), |1)}, the 
effective Hamiltonian becomes 



H (r — 



-E/2 
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where F(t) is the coefficient function in the exact O- 
operator D(t,z*) = F(t)L ( 26]) and satisfies 



d t F(t) 



7 



+ {- 1 + iE)F + F 2 



(5) 



and F(0) = 0. When the initial state is chosen as \i[>o) = 
(1/a/2)(]0) + |1)), by Eq. © and in the rotating frame 
of H sys , the fidelity measuring the survive probability of 
the initial state under the dissipative dynamics is easily 
obtained as 



l + K[e-f° dsF ^] 



(G) 



where 7Z[-] stands for the real part. 

For the two-qubit model, the exact O-operator con- 
tained in the stochastic Hamiltonian (flj can be explicitly 
derived ^ as <3(i, z*) = F 1 {t)0 1 +F 2 (t)0 2 +iU{t 1 z*)0 3l 
where O x — L, 2 = afa B + a^af , 3 = a^a B , and 
U(t, z*) = f Q dsll (t, s)z*. These coefficient functions sat- 
isfy 



d t F 1 (t) = j + (- 1 + iE)F 1 +Fl 



3F| 



:U, 



d t F 2 (t) = (-7 + iE)F 2 - Fl + ^F X F 2 + F 2 2 - -U, 
d t U(t) = -2ijF 2 + (-2 7 + 2iE)U + 4FiC7, 



(7) 



with the notation U(t) = f Q dsa(t, s)U(t, s). The bound- 
ary conditions are given by Pi(0) = F 2 (0) — U(0) = 
and U(t,t) = -AiF 2 (t). 

It is advantageous to choose {|s), \a), |11), |00)} as 
basis, where |s) = (l/\/2)(|10) + |01)) and \a) = 
(1/\^)(|10) - |01>). If P is chosen to be spanned by the 
first two basis vectors |s) and \a), the stochastic Hamil- 
tonian can be written as 



Hrff — 
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V2i(z* + U z ) 
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-E J 



(8) 



where f(t) = Fi(t) — F 2 (t). From Eq. ([7]), we can check 
that 

d t f{t) = 1 - + {- 1 + iE)f + 2f\ (9) 

and /(0) = 0. Suppose that the initial state is the pure 
state |Vo> = P. (0)|«) + Pa(fl)\a) + gii(0)|ll) + g O o|00) 
with |P S (0)| 2 + |P a (0)| 2 = 1 and ?11 (0) = <z O o(0) = 0. 
From Eq. ©, we have P s (t) = e~ 2 ds f^P s {0), P a (t) = 
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P„(0) and = P s (i)|s) +P a (*)|a> + 9oo(*)|00). The 
concurrence |34j | for each trajectory is given by 



F(co t=lO) 



C(^) = | e - 4 /o^/Wp2(0)-P o 2 (0)|. 



(10) 



In fact, Eqs. © and (TTTfl) do not contain stochastic 
noise explicitly. By Novikov theorem, ensemble aver- 
age M[x(f)] = x(f). Thus M[C(^t)] = Ctyt). The 
concurrence C(p t ) = \p ss - p Qa - Psa + p as \, where 
Pss = M[\P s (t)\% Paa = M[|P a (0)| 2 ], and p sa = p* as = 
M[P*(t)P a (0)]. Straightforward calculation shows that 
C(pt) = M[C(ipt)}- Moreover, C{pt) is independent on 
Poo and any coherence terms involving 1 00) . Thus for the 
initial states with qn(0) = 0, Eq. (jlOl) is valid and is 
always equal to C{pt). 

In a similar way, we get an explicit expression for 
the fidelity: T = (</>oH^o) = (V>o|M[|V t )<Vt|]|V'o) = 
M[(^t)(A\^o)}. Thus 



^(t) = |e- 2 J?^*W|pa( )| + |p! 



(11) 



It should be noted that if |P S (0)| = 1 and P o (0) = 0, 
i.e. |V>o) = |«), C(^t) = T{t) = e-4/o'rfsKI/Wl. In this 
situation, the fidelity can be used as a reliable measure 
for controlling entanglement. 

Non-perturbative control of fidelity and entanglement. 
- The exact subsystem stochastic master equation ([3]) 
is the starting point of our non-perturbative dynamical 
decoupling theory. The exact expressions in Eqs. ^ and 
([TO]) are determined by the differential equations given 
by Eq. ([5]) and Eq. respectively, and the solutions to 
those equations are dependent on two physical quanti- 
ties 7 and E. The parameter 1/7 is the correlation time 
of the environment and is non-tunable for a given envi- 
ronment. Our control is introduced through the system 
energy shift E(t) = w+e(t), where u is the bare frequency 
and c(t) is a time-dependent control. This is similar to 
coherent control theory, where phases of external field 
are the control parameters. In principle, the control c(t) 
can be an arbitrary function of time that can be realized 
by an external field applied to the system, including rect- 
angular pulse sequences with arbitrary strength and fre- 
quency. The zeroth-order perturbation or idealized pulse 
approximation shows that if the frequency of the pulse 
is sufficiently high, the information leakage rate of sub- 
space to the environment or the rest part of the system 
could be efficiently suppressed. In this Letter, we will go 
beyond perturbation and the idealized pulses to explore 
the non-perturbative regimes and control thresholds. 

The control parameters consist of the period of the 
rectangular pulse r, the duration A and strength ^> / A, 
respectively, i.e. c(i) = 'J'/ A for regions jit- A < t ^ nr, 
n ^ 1 integral, otherwise c(i) = 0. The pulse frequency 
is determined by the dimensionless proportion parameter 
t/A. The smaller r/A, the more frequent pulses are 
applied. 




FIG. 1. The fidelity at the moment ujt — 10 of the single- 
qubit system under control. It is initially prepared in |^o) = 
(l/\/2)(|l) + |0}). The bottom triangular area (dark blue) 
has no physical meaning since r is always larger than A by 
definition. We choose ^ — uj and 7 = 0.5. 



C(co t=10) 




FIG. 2. The concurrence at the moment ojt = 10 of a two- 
qubit system under control. It is initially prepared in l^o) = 
I a). Again, the bottom triangular area (dark blue) has no 
physical meaning since r is always larger than A by definition. 
We choose $ = uj and 7 = 0.5. 














t=2A 






t=3A 






t=4A 







FIG. 3. The concurrence dynamics of a two-qubit system 
prepared in |i/>o) = \s) with different A's and r's. (a) uiA = 
10" 2 ; (b) ujA = 2xl0~ 2 ; (c) ojt = 5xl0" 2 ; (d) cur = 6xl0" 2 . 
We choose $ = uj and 7 = 0.5. 
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FIG. 4. r/A — ^ phase diagram of the concurrence at the 
moment ujt = 10 of a two-qubit system under control. It is 
initially prepared in | i/^o) = I s )- We choose A = 10 -2 /u> and 
7 = 0.5. 
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FIG. 5. The concurrence damping time T(C = 0.95)/(wt) of 
entanglement of a two-qubit system prepared in |t/jo) = \s) vs. 
7 with different r/A's. We choose A = 0.02cj£ and >]/ — ui. 

In what follows, we will explore the parameter regions 
that allow efficient control of quantum coherence with 
more practical pulses. For the single qubit model, we plot 
the contour of the fidelity at uit — 10 as the function of 
the pulse period r and the pulse duration A, which may 
be called a r — A phase diagram. The result of ideal pulse 
sequence (A — > 0) is at the left-bottom corner in Fig. (p} . 
Surprisingly, our calculations show that there is a large 
parameter region for the non-idealized pulses, where the 
initial subspace or state \ipo) = (l/\/2)(|l) + |0)) can 
be protected with the same high fidelity as the idealized 
pulses, as shown in the red zone of Fig. ([!} and charac- 
terized byl<r/A<3. More transparently, we see that 
the fidelity remains above 0.95 at time uit = 10 for all 
the parameters satisfying 1 < r/A < 3. However, in the 
parameter region characterized by r > 3 A, the effective 
control becomes impossible. It turns out that what mat- 
ters for dynamical control is the ratio r/A rather than 
A or r alone. 

For the two-qubit model, our initial state is \ipo) = I s )- 



Without control, the state will lose its entanglement in 
time. Fig. ([2]) plots r — A phase diagram of concurrence 
at the time point uit = 10. The results clearly show 
there exists a threshold that dictates where the efficient 
entanglement control is possible (C 0.95). More specif- 
ically, the entanglement can be protected in the region 
of 1 < r/A < 8/3, we see again that the idealized-pulse 
approximation is only a small portion of the large per- 
missible parameter region. 

Fig. (j3|) shows the control condition of dynamics. In 
Fig. ©(a) and (b), we fix A and use different r values to 
show the effect of r; while in Fig. (0(c) and (d), we show 
the effect of A by fixing r. For a long time scale uit ^ 30, 
each value of A constrains the value of r to yielding a 
good entanglement control [that is C(t) is close to C(0)] 
and vice versa. We therefore conclude that as long as 
r/A is not too large, a quantum state or its entanglement 
can be well protected. This is the only condition for 
dynamical control to work. In addition, we have checked 
the pulse strength [See Fig. (jH)]. Interestingly, we see 
that the pulse strength does not play an essential role in 
quantum coherence control in contrast to the idealized 
pulse case where the the pulse strength is assumed to be 
effectively infinite. 

Note that Figs. (HJ-© are plotted with 7 = 0.5, which 
corresponds to a strongly non-Markovian regime. Fig. ([5]) 
is used to show how the environmental memory times 
affect our quantum control scheme. We use T to denote 
the time after which the concurrence decays below 0.95. 
Fig. ([5]) is plotted in a logarithmic scale to give rise to 
a better view of the entanglement change against 7. We 
emphasize that the non-Markovian bath is essential for 
the effective quantum control as shown in Fig. ([5]). 

Conclusion — We have established a non-perturbative 
approach to dynamical control theory based on the Fesh- 
bach PQ-partitioning technique and the non-Markovian 
quantum trajectory. Surprisingly, we have found the 
threshold and parameter regions of non-ideal pulse se- 
quence where dynamical decoupling can be used to pro- 
tect the entanglement and fidelity of the subsystem. We 
show that the idealized pulse approximation merely occu- 
pies a tiny part of the large permissible parameter region 
that is able to suppress the noise. This suggests that the 
zcroth-order approximation may not be necessary and 
thus prompts a more flexible experimental implementa- 
tion of the dynamical decoupling scheme. Additionally, 
we also testify to the environmental memory effect on the 
decoupling process. 
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